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1. Introduction 


Plasmopara viticola is the causal agent of the downy mildew, the most severe disease of the 
grapevine leading to economic damages (Wong et al., 2001). In order to prevent downy mildew, 
fungicide treatments are required, but they are dangerous for the environment and human health 
(Kab et al., 2017). Optimal scheduling and selection of treatments is the key to managing 
downy mildew in an eco-friendly way (Chen et al., 2020). This goal is quite difficult to achieve 
due to the variability shown by downy mildew among years. Indeed Plasmopara viticola growth 
mostly depends on variables like temperature and rain, plant’s genotype and soil conditions. The 
latter are usually assumed to be homogeneous in the considered vineyard, possibly because of 
the difficulty in obtaining local measurements. Meteorological variables are typically measured 
at whole-field levels, despite that Plasmopara viticola growth depends on microclimate (Bove 
et al., 2020a). Simulations of the key steps in the biological process of the pathogen have 
been performed to obtain information about airborne sporangia, sporangia availability, relative 
severity and number of lesions in secondary infection cycles (Brischetto, et al., 2021) (Bove et 
al., 2020b). Unfortunately these important deterministic models do not also provide information 
on the variability of the above attributes describing events related to the infection. 

In this work, we propose a Bayesian prior-predictive approach (Gelman, et al., 2017) where 
future environmental conditions and the probability of infection both depend on the selected 
treatment. A multi-attribute utility function taking the three most important variables as argu- 
ment has been elicited to describe the utility of consequences following the decision to treat 
the vineyard (Lavik, et al., 2020): the expected values under alternative decisions enable the 
winemaker to take the optimal decision of treating the vineyard or not. 


2. Methods 


In this section the approach followed to support the decision maker is described. 


2.1 Scenarios 


In this study intervals of temperature values and of humidity promoting the disease were defined 
by exploiting the information available in the literature. The following scenarios were defined: 
(i) a temperature favorable for pathogen’s growth but not for humidity, (Temperature > 10°C 
and < 30°C, Humidity < 0.8) labeled as ’Useful, N-Useful”; (ii) a temperature not favorable 
for pathogen’s growth and a favorable humidity (Temperature < 10°C or > 30°C Humidity 
> 0.8), labeled as ”N-Useful, Useful”; (iii) a temperature and humidity both favorable for 
pathogen’s growth, labeled as ’Useful, Useful’(Temperature > 10°C and < 30°C, Humidity 
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> 0.8); (iv) neither temperature nor humidity favorable for pathogen’s growth (Temperature 
< 10°C or > 30°C with Humidity < 0.8), labeled as *N-Useful, N-Useful’. Given that scenario 
ej (j € {1, 2,3, 4}) is realized in the vineyard, the expert must take the decision ”to treat”, a1, 
or not to treat”, ag 


2.2 States, actions, consequences 


Expected values of the probability 7; ; of infection for one leaf sampled from the vineyard given 
each environmental scenario e; and decision a;,i € {0, 1}, were elicited under the assumption 
that all of these combinations of temperature and humidity lasted from dawn to sunset just 
before taking the decision. After assuming that (m; ; | e;,a;) = Beta(a;,;, Bij), the values of 
model parameters a; ; and 6; į were defined for each pair scenario-treatment 7, j by fitting a Beta 
distribution to the elicited quantile 0.9 and the elicited expected value of 7;,; given aj, €j, i.e. 
pairs made by an action and a temperature-humidity scenario (Table 1). The implied credible 
intervals were checked by the expert (Table 1) without finding any need of refinement. 

Higher levels of variability characterize the prior-predictive distribution under no chemical 
treatment (ao) in comparison to the decision of treating (a;). In Table 1, the expected value 
of the probability of infection is shown for each scenario, p(741 | a;,e;), together with other 
elicited quantities. 


Table 1: Elicited expected values of the probability of infection in the considered scenar- 
os; ’Useful” (’N-Useful’”) means able (unable) to produce the infection; T=Temperature and 
H=Humidity. 


Treatments Scenarios e,,...,¢4 Probability Credibility Parameters 
{ao, a1} T H Elni5] Interval: 0.8 a, 5 Bij) 

0 Useful N-Useful 0.75 0.67296, 0.80032) 40.50, 13.50) 

N-Useful Useful 0.70 0.62413, 0.74968) 43.17, 18.50) 

N-Useful N-Useful 0.06 0.00066, 0.10263) 

) 
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( 
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Useful Useful 0.80 (0.72362, 0.84969 

Useful N-Useful 0.50 (0.46957, 0.52000 221.50, 221.50) 
( 169.33, 254.00) 
( 14.89, 134.00) 
( 


50, 112.50) 


N-Useful Useful 0.40 0.3696, 0.42000) 
N-Useful N-Useful 0.10 0.06991, 0.12001) 
Useful Useful 0.30 0.2696365, 0.32002) 
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Two attributes were defined to quantify the impact of a selected treatment on soil and bio- 
diversity of the vineyard at the subsequent time point t + 1 (e.g. next week) after the decision- 
action: 


e S41 : a score that classifies the degree of cleanness of soil after chemical treatment (in- 
cluding derived side products), Qs,,, € {1,2,...,5}, where s1 = 1 for the worst state 
after 10 years from treatment, and s;,, = 5 for the cleanest case after 10 years; 

e biy : a biodiversity score to classify the degree of biological diversity, (,,, € {1,2,...,5}, 
thus b,,; = 1 refers to the worst state of biological diversity after 10 years from treatment 
and b:4; = 5 is the best diversity class after 10 years from treatment. 


Given that the winemaker is willing to consider the two attributes on equal footing, a value 
function averaging and rescaling biodiversity and soil scores was considered as an environmen- 
tal summary of the future state: fc o2+1 = ((St41 + br41)/2 — 1)/4, with Qs = [0, 1]. In order 
to recognize the inherent uncertainty of f,4 141, a prior distribution was elicited by restricting 
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Figure 1: Contour plot of the utility function. 


the attention to the decision of treating, p( fs »t+1 | a1) = Beta(¢1, 2), because the decision 
of no treatment ao is associated with no change of biodiversity and nor of soil: a degenerate 
probability distribution follows under ag. For this reason the value of fs + was also calculated 
at the time of decision, thus p(fs,5,+1 | a0) = Tf, „+( f). The elicited value of the two parameters 
is d; = 57, 2 = 22, thus the treatment has a medium impact on the environment (quantile 0.1 
of fs» is 0.6559175; quantile 0.9 of fs b+ is 0.7846756). Hereafter, the probability of healthy 
leaves T; ; = 1 — 7;,; will be considered in the utility function. 
Under conditional independence of future attributes, the prior predictive distribution is 


P(fobt+1» Tij | Ís bt 1, G2, €j, @) = 
Beta(Tij | aij, Bij): [Betal fs vtri | $1, 92) Tla) + Irol f) Jo(a)] 0) 


thus the expected value of the utility function U (fs bt+1, Tij) is 


E|O (fs bt+1, Tij) | ai, ej] = ATES P(fs,bt+1, Tij | fsbt 1, Q2, €j, ai) dO 
6 


where @ is the vector of all model parameters. In the following, the current value of environ- 
mental summary is f,», = 1 under ap, i.e. a fully unmodified environment is in place. 


2.3 Elicitation of the utility function 


An utility function was elicited with arguments the environmental summary and the probability 
of healthy leaves: under mutually utility independence (French et al., 2000) (Keenye et al., 
1993): 


U (febir Tij) = kiU (fs bt+1) + k2Ua(Tij) +k ki ka Url fepe) U2lTiz) 


where k satisfies 1 + k = [J2 (1 + k ky); Ui(ai) = 9 Beta(z | pii, ~2«)dz,i = 1,2 are 
marginal utility functions which depend on parameters %1, ; and W,;; the best z¥ and worst q? 
cases take value equal to 1 and 0 respectively; the weights are elicited so that kı = u( fžp t41, T3) 
is the utility value associated to the best value for the environmental summary and the worst 
value for the probability of a healthy leaf; similarly, k2 = u(7j,, f24111) is the utility value as- 
sociated to the best value for the probability of a healthy leaf and the worst for the environmental 
summary. After eliciting U, and U> a graphical exploration was performed with the expert to 
check for the need of refinement (Figure 1). The optimal decision a! under condition e; follows 


from the expected values of the utility function: a = arg maxje4o,1; E[U(fs,oe41, Tij) | ai, ej]. 


235 


3. Results 


The expected values of the utility function were computed for each scenario as described in 
the previous section. In Table 2 the main results are shown. 

By comparing the different scenarios under different decisions, it was found that for e} = 
*Useful N-Useful’, the expected utility was higher in the ”not treat” case (a = 0), than ”treat” 
case; when e2 = ”N-Useful Useful’, the expected utility was higher in the ’treat” case (a = 1), 
than not treat” case; for e3 = ”N-Useful N-Useful’, the expected utility was higher in the ’not 
treat” case (a = 0), than treat” case; finally, when e4 = ” Useful Useful’, the expected utility 
was higher in the ’treat’” case (a = 1), than ”not treat” case. 


4. Discussion and conclusion 


Optimal scheduling and managing of treatments is a way to reduce the environmental impact 
of agriculture. This goal is quite challenging while dealing with phytopathogens that have high 
infectious potential and that may produce extensive and severe damage. Plasmopara viticola, 
the main enemy of viticulture, is one of these phytopathogens requiring the adoption of highly 
tuned prevention strategies. The wide adoption of treatments based on copper and sulphuric 
compounds is leading to over-accumulation in the soil, especially of copper, which causes a 
phytotoxic effect on the grapevine. They also have a negative impact on biodiversity by reducing 
the number of species and weakening the ecosystem in the long term. 

The optimal decision about treatment with chemicals rests on the available (prior) informa- 
tion about the risk of infection at decision time, the probability of observing a healthy leaf after 
treatment and the expected impact on the environment. The availability of data collected in the 
vineyard of interest is the natural next step to improve the performance of the decision process 
by better calibrating expectations and beliefs: here the advent of low cost sensors for oospores 
could lead to decisions taken for local microenvironments. Furthermore, agronomist’s prefer- 
ence scheme over prospects coded into the elicited utility function is crucial in order to define 
a trade-off between environmental sustainability and yield, both for quantity and quality. Here 
the four most fundamental scenarios of climatic conditions have been considered but a multi 
value discrete scale on more intervals for several other variables could increase the resolution 
of the description, when needed. Similarly, a direction for further research could be a more 
detailed description of both environmental changes and end products, grapes, by choosing key 
chemical components required to produce high valued wine. 


Table 2: Expected values of the utility function for each scenario considered; ’Useful” (N- 
Useful”) means able (unable) to produce the infection; T=Temperature and H=Humidity. 


Treatments Scenarios e1,...,e4 Expected Value of 
{ao, a1} T H Utility function 

0 Useful N-Useful 0.251 

0 N-Useful Useful 0.253 

0 N-Useful N-Useful 0.959 

0 Useful Useful 0.250 

1 Useful N-Useful 0.231 

1 N-Useful Useful 0.374 

1 N-Useful N-Useful 0.902 

1 Useful Useful 0.581 
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The proposed utility function was based on cumulated Beta distributions resembling to s- 
shaped curves. This is not the only possible choice, e.g. logistic functions could be used instead, 
as well as many other functions. Nevertheless, the fundamental feature that we believe should 
not change is the presence of high utility values only when high values are present both for 
the environmental attributes and for the leaves: this is quite expected in view of the increasing 
importance of environmental sustainability in agricultural decision-making processes. 

The end-user should not take the elicited functions as a black box reference ready to be 
exploited. The elicitation of soil and biodiversity classes is strongly dependent on the consid- 
ered vineyard and on the selected chemical, e.g. more or less impacting and more-less effective 
against Plasmopara viticola. Furthermore, our utility function could be extended to include 
more specific sustainability indexes, more attributes describing quality and yield of grapes, and 
even alternative types of chemical treatment. Any extension in the above directions should 
always put the individual preference scheme of the winegrower at the core of an unbiased elic- 
itation procedure. 
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